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We analyze the geometric structure and mechanical stability of a complete set of isostatic and 
hyperstatic sphere packings obtained via exact enumeration. The number of nonisomorphic isostatic 
packings grows exponentially with the number of spheres N ^ and their diversity of structure and 
symmetry increases with increasing N and decreases with increasing hyperstaticity H = Nc — Njso , 
where Nc is the number of pair contacts and Njso = 3 A/" — 6. Maximally contacting packings are 
in general neither the densest nor the most symmetric. Analyses of local structure show that the 
fraction / of nuclei with order compatible with the bulk (RHCP) crystal decreases sharply with 
increasing N due to a high propensity for stacking faults, 5- and near- 5-fold symmetric structures, 
and other motifs that preclude RHCP order. While / increases with increasing H, a significant 
fraction of hyperstatic nuclei for N as small as 11 retain non-RHCP structure. Classical theories 
of nucleation that consider only spherical nuclei, or only nuclei with the same ordering as the bulk 
crystal, cannot capture such effects. Our results provide an explanation for the failure of classical 
nucleation theory for hard-sphere systems of A/" ^ 10 particles; we argue that in this size regime, 
it is essential to consider nuclei of unconstrained geometry. Our results are also applicable to 
understanding kinetic arrest and jamming in systems that interact via hard-core-like repulsive and 
short-ranged attractive interactions. 
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I. INTRODUCTION 

Crystallization of monodisperse hard spheres is a com- 
plex problem for several reasons. In the absence of attrac- 
tive interactions, crystallization is associated with min- 
imization of free volume. For bulk systems, the ground 
states, the FCC and HCP lattices, as well as other stack- 
ing variants of hexagonal planes, possess volume fraction 
0xtai = tt/a/TS .7405 [T. In addition, there are an ex- 
ponential number of rigid packings with volume fractions 
that range from random close packing [2^ to ^xtai- The 
large number of metastable structures and large barriers 
separating the amorphous and crystalline states leads to 
formation of amorphous structures if the quench 

rate is not sufficiently slow. This makes understanding 
crystal nucleation in these systems particularly impor- 
tant since jamming or glass-formation may be avoided 
only through nucleation and growth of crystallites. Since 
the bulk crystal state maximizes ^, classical nucleation 
theory suggests that the densest packings of N spheres 
within a (minimal) spherical volume V are optimal nu- 
clei. However, this approach should work if and only if 
these packings possess the same structural order as the 
bulk crystalline phase. A recent study by Hopkins et al. 
[8] showed that this condition fails; the densest packings 
did not in general have FCC, HCP, or Barlow H [10] 
order. Instead, their surface order was dominated by 
the spherical boundary conditions, and so they may not 
correspond to the stable nuclei that form in unconfined 
geometries, i.e. within an arbitrary volume in a larger 
system. 

It is important to note that small nuclei may be dis- 



tinctly aspherical^ with a wide variety of shapes, symme- 
tries, and formation probabilites. Specifically, A^-sphere 
nuclei with Nc contacts can form A4{N^Nc) packings of 
distinguishable shape, symmetry, and entropy [TT]. How- 
ever, the full range of shapes, symmetries and statistical- 
gometrical properties of such packings has not been quan- 
titatively characterized, even for small N. There have 
been numerous studies of the phase diagrams and crys- 
tal nucleation and growth in systems of hard spheres and 
sticky hard spheres p!2HT7] . However, studies of hard 
sphere crystal nucleation in particular have shown that 
quantities inferred from experimental results and classi- 
cal nucleation theory can differ from simulation results 
[18 by orders of magnitude. Thus, there is a need to 
characterize the structural properties of nuclei posessing 
arbitrary geometry and interactions with the surround- 
ing fluid [17, 19 to gain a more quantitatively accurate 
description of crystal nucleation and growth. In this pa- 
per, we analyze the statistical-geometrical properties of 
small sticky-hard-sphere nuclei to lay a foundation for 
a more quantitative understanding of crystal nucleation 
in systems such as colloids with hard-core-like repulsions 
and short-range attractions [TTl [20] . 

Determining these properties is largely an exercise in 
the geometry of finite sphere packings ^2TH25] . The two 
key mathematical problems are "What is the maximum 
number of contacts N^'^^{N) that can be formed by N 
monodisperse spheres?" and "How many different ways 
A4{N^Nc) can N spheres form Nc contacts?" Solving 
these problems simultaneously yields complete sets of 
isocontacting packings of hard spheres and of isoener- 
getic states of sticky hard spheres. These in turn have 
applications to physical problems ranging from crystal- 
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lization and jamming |3H3 [T^^ TF, TH^ to cluster physics 
[IIl|22l|26^ to liquid structure [16 , 21-22, to protein fold- 
ing [30l [31] , as well as engineering applications such as 
circuit design [32] and error-correcting codes [33] , 

Despite this wide applicability, progress in obtaining 
solutions has been slow. Determining N^'^^{N) corre- 
sponds to the generalized Erdos unit distance problem 
[34 in three dimensions, which remains unsolved, while 
determining A4{N^Nc) has been proven to be algorith- 
mically "NP-complete" [35 . The latter condition im- 
pedes calculation of J\4{N,Nc) via Monte Carlo or re- 
lated methods [36]. Consequentially, M{N, Np''''{N)) 
has been determined for N as large as 10 only recently 
[22| l25l l37] . Here we present an efficient method for 
finding N^^^{N)^ M{N^Nc) and the permutational en- 
tropies of distinguishable sphere packings. A4 is an in- 
teger for isostatic {Nc = 3A^ — 6 = Njso) and hyper- 
static {Nc > 3N — 6) clusters of sticky spheres [21] [22] 
wherein each sphere possesses at least 3 contacts. We fo- 
cus on packings satisfying these necessary [38 conditions 
for mechanical stability since they correspond to solidlike 
clusters that likely play an important role in nucleation 
and growth of crystals. 

We find all isostatic and hyperstatic sticky hard sphere 
packings for < 11 via exact enumeration, and present 
novel analyses of several statistical-geometrical proper- 
ties of this complete set of packings that are relevant for 
understanding crystallization and jamming. Several dra- 
matic features are associated with the increasing maxi- 
mum hyperstaticity Hmax{N) = N^''''{N) - Njso for 
> 9. Key amongst these are that maximally con- 
tacting packings are in general quite different from the 
densest packings. Minimal energy (maximally contact- 
ing) packings are not necessarily either the most compact 
or symmetric. Instead, the most symmetric and compact 
packings are often mechanically stable "excited states" 
with Nc < N^^'^iN). While many of these stable pack- 
ings correspond to "on-pathway" nuclei possessing struc- 
tural order consistent with the bulk (Barlow-ordered) 
crystalline phase, many do not. The latter correspond 
to 'off-pathway" nuclei possessing structural motifs in- 
compatible with Barlow order, such as 5-fold-symmetries, 
stacking faults, and twin defects, all of which are known 
to impede crystallization [12i4T5l [2Q]. The fraction of iso- 
static nuclei possessing non-Barlow order grows rapidly 
with increasing N to nearly 95% for A/" = 11. Crystalline 
order increases with increasing hyperstaticity, yet a sig- 
nificant fraction (~ 50%) of hyperstatic nuclei for N as 
small as 11 retain non-Barlow structure. 

The outline of the remainder of the manuscript is as 
follows. In Section [TT] we describe our exact enumera- 
tion procedure, focusing particularly on advances beyond 
those employed in previous studies [22 , 25 , 37 . Section 
III presents analyses of the structure and symmetry of the 
sphere packings, including size, shape, and the statisti- 
cal prominence of key structural motifs. Section IV re- 



lates our work to other studies of crystal nucleation that 
include interactions of nuclei with the bulk fiuid and a 



range of nucleation pathways. In Section |V| we place our 
results in context with other recent work, and conclude. 
Finally, Appendix |A] describes mechanical stability anal- 
yses of the nuclei. Appendix [B] shows convergence of the 
exact enumeration procedure, and Appendix [C] explains 
implicit contact graphs for sphere packings. 



II. METHODS 

A key feature of (sticky) hard sphere packings is that 
their (isoenergetic) isocontacting states are in general 
highly degenerate. Many distinguishable arrangements 
( "macrostates" ) of N particles with Nc contacts are pos- 
sible. Here we employ a particular definition of the term 
macrostate. An A^-particle, A^c-contact macrostate is de- 
fined by a unique set of N{N — l)/2 squared interpar- 
ticle distances {r'^j} [38l|42]. Different macrostates have 
different {rfj} and "shapes". In general, systems of A^ 
spheres with Nc contacts possess Ai{N^Nc) distinguish- 
able macrostates. For example, A1(6, 12) = 2 since ex- 
actly two six-particle macrostates exist for systems with 
Nc = Ar^^^^(6) = 12 (Fig.[l} [21]. 

(a) ^ ^ (b) 



FIG. 1: (Color online) Macrostates: A^(6, 12) = 2. The oc- 
tahedral structure (a) has high symmetry and low permuta- 
tional entropy {uja = 15), while the capped trigonal bipyramid 
structure (b) has low symmetry and high permutational en- 
tropy (cjb = 180) [TT . The position of the green (rightmost) 
sphere in (b) implies a stacking fault. Note that in this and 
many subsequent figures, sphere sizes are reduced for visual 
clarity, and the connecting bars indicate pair contacts. 

There are many ways to organize indistinguishable 
spheres into any given macrostate; these correspond to 
permutations of particle indices {i^j} that preserve {rfj}. 
We refer to the number of allowed permutations as the 
number of microstates uj^ corresponding to a particu- 
lar macrostate k. Differing structure and symmetry of 
macrostates imply they have different Uk {i.e. permu- 
tational entropies). For example, the highly symmetric 
octahedral structure shown in Fig. [TJa) has uOa = 15, 
while the less symmetric structure shown in panel (b) 
has 002 = 180. Note that (a) is a subset of the FCC 
and HCP lattices, while the capped trigonal bipyramid 
(b) is a stack-faulted structure; for example, if the green 
sphere is removed, the remaining 5 spheres have HCP 
order. Such effects have important implications for nu- 
cleation, as macrostates with higher uj will form with 
greater probability Jl] 03] . These will be discussed in 
detail below. 
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The potential for sticky hard spheres with diameter D 
and contact attractions is 40 : 



Generate adjacency matrix 
from binary number 



Apply Geometric 
Rejection Rules 



Solve 
for{r} 



Analyze structure 
and stability 



Ussir) = 




(1) 



Since f/ss(r) possesses no scale, and we consider monodis- 
perse systems, the value of D is arbitrarily set to unity 
below. Hard-sphere constraints imply that the center- 
to-center distances Vij between unit spheres i and j with 
positions fi and fj obey Vij > 1, where the equality holds 
for contacting pairs. 

In this manuscript, we enumerate the global and low- 
lying local potential energy minima ^ of Eq. [l] It IS 
worth noting that at finite temperature T and zero pres- 
sure, systems interacting via the sticky hard-sphere po- 
tential (Eq. [T]) will have no persistent contacts since the 
range of attractive interactions is exactly zero. For this 
reason Baxter ^41^ introduced the adhesive hard sphere 
(AHS) interaction potential 



oo 



UAHs{r) 







r <D 
, D<r<rc (2) 
r >rc 



where r is a temperature-like parameter and Tc is the in- 
teraction range. The packings reported in this work have 
configurations ({r}) that are identical to the correspond- 
ing AHS packings in the r ^ and Tc ^ limits. While 
AHS systems have also been shown to possess thermo- 
dynamic anomalies in the Tc ^ limit [27 , several the- 
oretical studies [71 [TTl have shown that these vanish 
when Tc is as small as ^ .01 e.g by identifying Eq. 
[2] as the short-range limit of the attractive square-well 
potential. Applicability of our studies to systems inter- 
acting via short (but finite) range potentials is discussed 
in Section [mCl 

Employing an infinitely narrow potential well also al- 
lows all isostatic and hyperstatic A/'-sphere configura- 
tions to be conveniently characterized by N x N adja- 
cency matrices A with Aij = 1 for contacting particles 
and Aij = otherwise. In general, the configuration 
{r} = {ri,r2, ...r/v} can be solved (c/. Section IID) from 
A [38^ provided two minimal conditions for mechanical 
stability are met [44]: (i) each particle possesses at least 
3 contacts and (ii) Nc > Njso- Throughout this pa- 
per, we refer to these as conditions (i) and (ii). Using 
an efficient exact enumeration algorithm schematically 
depicted in Fig. [2j we identify all A and {r} that corre- 
spond to nonoverlapping configurations with Nc > Njso 
contacts. These configurations all have equal potential 
energy U = —^^j^iAij = —Nee for sticky spheres, and 
their {r} are identical to corresponding A/'c-contact hard 
sphere packings (to see this, consider approaching the 
limit e ^ from e < 0). Overlapping configurations with 
one or more Vij < 1 have infinite U . We refer to configu- 
rations with no interparticle overlaps as "valid packings" 
and overlapping configurations as "invalid packings." 



FIG. 2: Schematic of our exact enumeration algorithm includ- 
ing structural and stability analyses. 



The rest of this section describes our exact enumera- 
tion method in detail, following the scheme depicted in 
Fig.|2j For each N and Nc-, we perform complete enumer- 
ation by efficiently iterating over all adjacency matrices 
(Section |ll A ). Valid packings are found by applying geo- 
met rical and graph-theoretic rejectio n rule s (Sections |II B 
and |II Cj ) and then solving (Section II D ) for the struc- 
ture {F} of nonisomorphic packings passing these rejec- 
tion rules. The structure and stability of valid packings 
are then analyzed as described in Sections [II E| and [II F[ 



A. Exact Enumeration Method 

Since all elements of adjacency matrices are or 1, they 
correspond to binary numbers B. The matrices are sym- 
metric, and diagonal elements are zero by convention, so 
any A may be uniquely associated with one 
distinct B. All sticky sphere packings may be found by 
iterating sequentially [45 over the B and mapping each 
to an adjacency matrix . 

The number of adjacency matrices that must be iter- 
ated over to ffnd all macrostates and microstates for ffxed 
N and Nc but with no constraints on the arrangement of 
the elements {i.e. arbitrary topology) is 



arb 



[{N^-N)/2]\ 
Ne\[{N^ - N)/2- N, 



cj! 



(3) 



JVarb grows faster than exp(A^) and rapidly becomes pro- 
hibitively large; for example, J^arb = 3.824 • 10^^ for 
N = 11 and Nc = 27. 

However, all macrostates can be found with greatly re- 
duced computational effort through appropriate selection 
of "topological" constraints on the elements of A. Biedl 
et al. [46 proved that all connected sphere packings ad- 
mit linear polymeric paths, i.e. for any valid packing, one 
can always permute particle indices so that the packing is 
fully traversed by a "polymeric" A with A^^^+i = 1 for all 
i. As in Ref. [37], we impose polymer topology by fixing 
= 1- Thus — 1 elements of A are fixed to unity, 
and the remaining N{N — 1) /2 — {N — 1) elements are left 
unconstrained. This arrangement reduces the number of 
binary numbers and adjacency matrices over which one 
must iterate to 

^ [{N^-^N^2)/2]\{\Nc-{N-m-' 

[{N^ - 3^- + 2)/2 - {{Nc -{N - 1))]! * ^ ^ 

For each N and A'c, we iterate sequentially over the Mpoi 
binary numbers and adjacency matrices as illustrated in 
Fig. [4j Mpoi/Marb decreases faster than exponentially 
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with increasing with a corresponding reduction of 
computational effort. For the purposes of calculating 
M.{N^Nc)^ this provides a speedup of about 3 orders of 
magnitude for the largest TV considered here (Fig. |3]) 




10 11 



FIG. 3: (Color online) Computational effort reduction ob- 
tained by using polymeric topology in enumeration. The 
graph shows faster than exponential decay of Mpoi / Marh with 
increasing N for Nc — 3N — 6. 

Here we are considering colloidal clusters with no fixed 
topology. Enumerating over adjacency matrices with 
polymer topology naturally produces {cok} correspond- 
ing to packings with polymer topology, i.e. different 
absolute macrostate populations uJk and ratios oJk/^j 
{j^k G {1,A^}) due to entropic factors such as block- 
ing [37]. However, the uj^ for colloidal clusters may be 
calculated via symmetry operations: 



CkN\/Ak 



(5) 



where Ak is the number of automorphisms of the adja- 
cency matrix corresponding to macrostate k [47 . C/e = 2 
for macrostates possessing chiral enantiomers and 1 for 
those which do not |^. The total number of mi- 
crostates for nuclei with N particles and Nc contacts is 
then 



MiN,N,) 



(6) 



k=l 



Note that while Eq. [s] treats particles {l^N} as distin- 
guishable [48 , we have verified that Eqs. [5][6] produce 
the same uj^ and Q produced by an alternative method 
pT| [25] that treats particles as indistinguishable and cal- 
culates UJk using the number of symmetries c possessed 
by the coordinate solutions {vk}- 



B. Geometric Rejection Rules 

Valid sticky sphere packings correspond to A^-vertex, 
A/'c-contact unit distance graphs that are embeddable 
[35] in three dimensions. A key advance in determin- 
ing embeddability of small packings was recently made 
by Arkus, Brenner and Manoharan [22l[25]. They used 
concepts from sphere geometry to develop geometric re- 
jection rules identifying invalid packings with patterns 



Generate adj. mat. 
for binary number B 



Apply M-level rejection 
rules to all M-subgraphs 
of adj. mat. 



M^M+l 




Store adjacency matrix 



NO 



B<Npo|(N,Nc)? 



YES 



Finished enumeration: identify 
nonisomorphic adj-mats 



FIG. 4: Schematic of Part 1 of our enumeration protocol: 
Geometric rejection rules. The sequential enumeration over 
binary numbers and application of rejection rules correspond 
to the left two boxes in Fig. [2] and are implemented as de- 
scribed in Sections III Alill CI 



within adjacency matrices. Geometric rejection rules fa- 
cilitate connections to graph theory and enable formula- 
tion of rules in terms of Boolean satisfiability conditions. 
These conditions can be conveniently organized into a 
series of "M" -rules that reject an invalid A based on pat- 
terns within M X M subgraphs of A, or restrict the way 
additional spheres can be added to form valid (M + 1)- 
particle packings. We apply these in order of increasing 
M to reject invalid packings. 

As illustrated in Fig. [4| we apply the rejection rules to 

all 



F{M) 



M\{N -M)\ 



(7) 



M xM subgraphs of each A. Any subgraph violating any 
rule indicates an invalid packing. For example, at each 
rule- level M, we test for over-connected "O" -clusters as 
follows. 0-clusters are defined by the M-subgraph and 
the O — M particles that contact at least one parti- 
cle in the M-subgraph. We independently verfied that 
^max(Q) is the same as found in Refs. [Ulllll for O < 11. 
If the total number of contacts in the O-cluster is greater 
than N^'^'^iO) {i.e. 30 - 6 for O < 9, 30 - 5 for O = 10, 
and 30 — 4 for O = 11), we reject the adjacency matrix. 
These overconnected-subcluster rules eliminate many in- 
valid packings and improve the efficiency of the code. 
Refs. ^22| [25] reported a complete set of rejection rules 
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for packings of N < 7 spheres. This is equivalent ly a com- 
plete set of rejection rules for M < 7 subclusters (with 
Mc = 3M — 6 contacts) within larger packings. We have 
extended this set to reject all invalid isostatic and hy- 
perstatic packings of < 9 particles, and many invalid 
packings of N > 9 particles. Note that in contrast to Ref. 
[25j, we do not use the triangular bipyramid rule [50 , 
nor explicitly search for conflicts in the distance matrix 
D {Dij = Tij) arising from different M < N subgraphs. 
Instead we employ geometric and graph- theoretic rejec- 
tion rules that do not require calculation of unknown dis- 
tances. Space constraints preclude describing our com- 
plete set of rejection rules for M < 9; here we highlight 
several rules not contained in Refs. p2 | |25]. 



Several rejection rules are obtained from known graph- 
theoretic results for the embeddability of sphere pack- 
ings. Kuratowski graphs [51 Km,n have m -\- n vertices 
and mn edges, with each of the vertices having degree n. 
The graphs i^3,3, i^3,4, and K^^4 are not embeddable as 
3D sphere packings [52l|53 ; these results imply rejection 
rules for (M = 6, Mc = 9), (M = 7, = 14), and 
(M = 9, Mc = 18) subgraphs, respectively. 



We make use of the fact that no structures with BCC 
symmetry are among the isostatic or hyperstatic pack- 
ings for sticky hard spheres [54 since placing a ninth 
sphere inside a cube implies overlap. Figure [sja) shows 
a cubic structure with M = 8 and Mc = 12. Placing a 
ninth sphere in the interior of the cube to form a putative 
A^ = 9, A^c = 20BCC packing implies an overlap of at 
least {V3 - l)/4, i.e. an < 1 - (Vs - l)/4). As illus- 
trated in Fig.[5j many M = 8 rejection rules are obtained 
by observing that a ninth sphere must lie in the interior 
of a cube or sheared cube if it contacts more than 4 of 
the 8 particles; any such placement implies particle over- 
lap. These rules eliminate many invalid N = 9, A/'c = 21 
packings not eliminated by M < 8 rules, and become 
increasingly effective at eliminating invalid packings for 
N>9. 



Other M = 8 rejection rules relate to "irregular" 
seeds lacking any underlying cube or sheared-cube topol- 
ogy. Figure |6ja) shows an M = 8, Mc = 17 packing 
that cannot be "4-kissed" to form a 9/21 packing. A 
ninth monomer cannot contact the four (blue and green) 
monomers because doing so would imply overlap be- 
tween the green monomers {i.e. an interparticle distance 

tg 0.615), indicated by the dashed red line. Figure 
b) shows an M = 8, Mc = 13 packing with the topol- 
ogy of a partial icosohedron. A ninth monomer cannot 
contact all eight to form a N = 9^ Nc = 21 packing with- 
out implying overlap. The rules shown in Figs. [5][6j to- 
gether with a few additional similar rules for M < 8, are 
sufficient to reject all invalid A/" = 9, A^c = 21 packings. 



C. Systematic development of additional rejection 

rules 

The number Q of rejection rules required to reject all 
invalid A^-sphere packings is expected to grow exponen- 
tially with N [50l |55]. Development of geometric re- 
jection rules "by hand", as described in the above sub- 
section, therefore becomes increasingly difficult as A^ in- 
creases. Here we report a systematic method for develop- 
ing additional rejection rules. We employ a "deep-seed" 
elimination procedure: 

1) Find all nonisomorphic "M-seed" graphs of M ver- 
tices and Mc < 3M — 6 edges satisfying minimal rigidity 
condition (ii) and passing all L < M rejection rules. 

2) Determine which of these can form known- valid 
packings of P = M+1 particles and Pc contacts by exam- 
ining all possible arrangements wherein the P^^ sphere 
contacts (Pc — Mc) spheres of the M-seed. 

3) A seed graph that can never have such an (Pc — Mc) 
kisser often yields a novel rejection rule. 

For example, we find 540 nonisomorphic M = 9, 
Mc = 20 seed graphs that satisfy condition (ii) and pass 
all M < 8 rejection rules. A tenth sphere can contact four 
(of the nine) particles to form a P = 10, Pc = 24 packing 
for only 197 of these seeds. The remaining 343 seeds can- 
not be a subgraph of any valid 10/24 packing. Therefore 
all A containing subgraphs isomorphic to any of these 
343 and a 10th particle contacting four of the particles 
in A^ correspond to invalid packings and are rejected at 
the M = 9 level. This "no-4-kisser" rule is particularly 
effective, eliminating 132 nonisomorphic invalid A^ = 10, 
A^c = 24 packings that passed all previously implemented 
rules. Most of these eliminated packings are invalid be- 
cause the putative 10th sphere contacts an open 4-ring on 
the surface of a 9-sphere seed. The additional four con- 
tacts can be formed without producing overlap only if the 
planar angle satisfies 27r/3 < < tt (Fig.[7|). While all 
M = 9, Mc = 20 packings fail to satisfy condition (i), 
many are insufficiently floppy for to fall within this 
range, i.e. addition of a 10th sphere implies < 27r/3 or 



D. Efficient Euclidean structure solver 

We solve for the Euclidean structure {r} of noniso- 
morphic packings which are not eliminated by any of the 
geometric rejection rules. The adjacency matrix yields a 
set of A'c equations and A^(A^ — l)/2 — A^c inequalities for 
ie[l,N],j e[i + l,N]: 

|fi-r,f =1 ; Aij = l 

(8) 

\n-rj\^ > 1 ; Aij = o. 

Solutions to Eq.|8]are valid N/Nc packings except in the 
case where they possess "implicit" contacts correspond- 
ing to the "=" case of the ">". Proper accounting of 
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^^^^ ^^^^ ^^^^ ^^^^ 

FIG. 5: (Color online) Schematic for M = 8 cube and sheared cube-based geometric rejection rules. No 9th particle may 
contact more than four particles of any subgraph isomorphic to the cube (panel (a)), sheared cubes with contacts across one 
(panel (b)), two (panel (c)), or three (panel (d)) faces. Similar exclusions apply to other subgraphs similar to panels (c-d) but 
with different topology, e.g. two cross-face contacts on opposite as opposed to adjacent faces. Panel (e): no 9th particle may 
contact more than four particles of the square antiprism with Mc — 16 or the sheared square antiprism with Mc — 17 (the 17th 
contact is indicated by the green line). 




FIG. 6: (Color online) "Irregular" M = 8 rules. Panel (a): 
Contact of a ninth monomer with the red and blue (dark- 
est shaded) monomers implies overlap of the green (lightest 
shaded) monomers (red dashed line, dgg < 1) and/or a planar 
angle ^ < 27r/3 (see Fig. |7|. Panel (b): Contact of a ninth 
monomer with each of the eight shown implies at least one 
overlap. 




FIG. 7: (Color online) A useful rejection rule for M-particle 
packings containing open 4-rings on their surfaces is that an 
(M+1)** particle can contact each of the 4 particles in the ring 
only if the planar angle (indicated by the arrows) satisfies 
27r/3 <ip <7T. 



attempts to find the roots of 

m = Ej>i S (Aij - 1) (|rl - r-f - l) + 
REj>i 5 {A.j) {\n -f/-l)e{l-\r-l-fj\)= 0, 

(9) 

where F{{f}) is the "error" function, 5 is the Kronecker 
delta function, and 6 is the Heaviside step function with 
0(x) = for X < and 1 for x > 0. The first term in 
Eq. [9] enforces non-overlapping contact between particle 
pairs with Aij = 1 and the second term is a repulsive 
term penalizing overlaps for particles with Aij = 0. For 
R = 10 the combination of repulsive force and step size 
control gives a large (order-of-magnitude) speedup over 
a version lacking these features. 

The iterative nature of our solver is illustrated in Fig. 
|8] Solutions are considered converged and a valid pack- 
ing is found when \F{{r})\ < TOL. If \F{{r})\ > TOL 
after nj Newton iterations, the solution is discarded and 
the process begins with a new {vinit}- If a solution is 
not found after ut attempts, A is rejected as an invalid 
packing. We find that failure of convergence of the struc- 
ture solver to converge to |F({r})| < TOL within ut at- 
tempts is sufficient to reject invalid packings, provided ut 
is sufficiently large. This sampling over different {vinit} 
is an important part of our enumeration procedure since 
the set of geometric rejection rules remains incomplete 
for > 9; further details are given in Appendix [B] 



E. Dynamical matrix analyses 



implicit contact graphs is key to exact enumeration stud- 
ies, both for determining M and for developing graph- 
theoretic rejection rules; see Appendix [C] for a discussion 
of these issues. 

We solve Eq. [S] efficiently using a multidimensional 
Newton solver with step size control ^ schematically 
depicted in Fig.jsj Initial conditions for the solver {rinit} 
are generated by placing N particles randomly within a 
cube of length A/", centered at the origin. The solver then 



Determining mechanical stability of packings is of great 
interest since packings with "floppy" modes have high vi- 
brational entropy. We determine stability using dynam- 
ical matrix analyses [44]. The Hessian matrix 



_dHJ_ 

dfidfj 



(10) 



has 3A^ — 6 positive eigenvalues for mechanically stable 
packings, but fewer for floppy packings. Since Eq. [l] is 
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to next adj. nnat. 



FIG. 8: Schematic of Part 2 of our enumeration protocol; 
Euclidean structure solver. This diagram corresponds to the 
third box in Fig. |2] All symbols (nr, nj, F{x), TOL) are 
described in the text. 



singular at r = 1, we (following Ref. [37 ) replace it by 



Uharm ('^) 



-e+ y (r - 1) , r < Tc 



(11) 







with U = ^jy^Uharmi'f^ij)' Note that Eg. [TTj reduces 
exactly to the sticky sphere potential (Eq. [T]) in the limit 
kc oo. We choose kc = lO^e and Tc/D = 1 ^/2e/k. 
For this kc, in all cases (since implicit-contact graphs are 
eliminated; see Appendix [C|), only the Nc pairs specified 
by A interact via Uharm- 



compatibility with the bulk crystal, and the presence of 
various structural motifs within nuclei. 

Point groups provide a convenient means of classifying 
sticky sphere packings. Macrostates with higher sym- 
metry have lower permutational entropy [11 and are of- 
ten associated with higher crystalline order. We evaluate 
point group symmetries of packings using the Euclidean 
solutions for {r} and the symmetry e valuator of Lee and 
Shattuck [57]. 

Barlow packings are hard-sphere packings com- 
posed of layered hexagonal-close-packed planes; their 
three-dimensional order may be FCC, HCP, or mixed 
FCC/HCP, but they possess no defects {e.g. stack 
faults.) They are optimal nuclei for hard- and sticky- 
hard-sphere crystals since they possess the same order- 
ing as the bulk equilibrium crystals, and have (j) 
tt/a/TS in the N ^ oo limit. We identify nu- 
clei with Barlow order by verifying that all {r'fj} are 
equal to values found in Barlow packings, i.e. {r^j} ^ 
{0, 1, 2, 8/3, 3, 11/3, 4, 5, 17/3, 6, 19/3, 20/3, 7, 22/3, 8, 
25/3, 9, 29/3, 10, 31/3, 32/3, 11, 34/3, ...} for all i and j. 

Structural motifs relevant to nucleation, e.g. 5- and 
near-5-fold symmetric structures and stacking faults, are 
identified through the presen ce of M x M subgraphs A' 
uniquely associated (Section II D) with the correspond- 
ing structures within the N x N adjacency matrices. For 
motifs associated with a pattern X, we identify the num- 
ber of macrostates A4x and fraction of microstates fx 
including these patterns via 



Mx{N,Nc) 



MiN,Nc) 



(12) 



where Gk{X) is 1 if structure of the k^^ macrostate 
matches the pattern and otherwise, and 



M{N,Na) 
k=l 



(13) 



where {uj} and Q are given by Eqs. 5][6 



We choose to identify these motifs as described above 
rather than alternatives such as determining the number 
of "crystal-like" particles as is common practice in the lit- 
erature [58] . The latter practice is better suited to studies 
of bulk crystallization, whereas we consider small nuclei 
where surface effects dominate. For example, while many 
studies have examined formation of nuclei with Barlow 
(FCC, HCP, and RHCP) ordering, we are not aware of 
any previous studies that quantitatively examined the 
fraction f Barlow of uuclci possessing such order. 



III. RESULTS 

F. Structural Analyses 

We now report results for the number, structure, sym- 
We analyze the structural order of nuclei using sev- metry, and dominant structural motifs within all iso- 
eral measures: crystallographic point group symmetry, static and hyperstatic sticky sphere packings for A/" < 11. 
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Results for M{N^Nc) are shown in Table |lj Values for 
< 10 are the same as reported in Refs. [25] [37] . 
A4 grows exponentially with increasing N. Exponential 
growth is expected for systems with short-range inter- 
actions and liquidlike order [59 , but has not previously 
been conclusively demonstrated for sphere packings [60] . 
The arguments by Stillinger and Weber in Ref. [59^ ap- 
ply to "large" N; our results suggest that sticky sphere 
packings are already in the large- A/' limit for > 9. 

To examine the degree of crystalline versus liquidlike 
order, we report the number of macrostates with Ci point 
group ordering and Barlow ordering {Mci and Msariow 
respectively). We also report the fraction of microstates 
fci and f Barlow with Ci and Barlow order. For isostatic 
packings, fd increases rapidly with while f Barlow de- 
creases rapidly for N > 7. Fractions of nuclei with Bar- 
low order increase sharply with hyperstaticity, consistent 
with the onset [4 of crystallization for Nc > Njso- 

The arguments of Phillips and Thorpe that glass- 
formation is optimized when systems are isostatic have 
been supported by many studies, including studies of sys- 
tems interacting via central forces. Our results support 
these arguments. Q decreases sharply with increasing hy- 
perstaticity (e.g. n{N = 11, Nc = 29)/n{N = 11, Nc = 
27) = 8.56 • 10~^), indicating an entropic barrier to in- 
creasing Nc beyond Njso- A large fraction of isostatic 
nuclei have liquid-like symmetry yet are solidlike in char- 
acter {i.e. mechanically stable [38 -see Appendix[A|) , and 
cannot change structure without breaking bonds, indicat- 
ing (for sticky spheres) an enegetic barrier to increasing 
Nc beyond Njso- These results provide a quantitative (if 
partial) explanation for earlier reports of glass formation 
by kinetic arrest in sticky hard sphere systems [71 |20l [62] . 
Specifically, both energetic and entropic barriers should 
impede nucleation of more ordered crystallites; our quan- 
tification of these effects may help explain why classical 
nucleation theory breaks down for ^ 10 [18 . In the 
following subsections, we will examine the shapes, sym- 
metries, and relevant structural motifs of sphere packings 
with A^ < 11 in quantitative detail. We found 99% of 
these packings to be mechanically stable; results of our 
stability analyses are discussed in Appendix [A| [63] . 



A. "Bulk" measures of shape and symmetry 



The complete set of packings reported here | 
hibits a great diversity of symmetries and shapes. 
[9|a) shows values of i?^. 



1 ^ 

N ^ 



N 



(r-)r 



^ ex- 
Figure 



(14) 



for all 9 < A/" < 11 macrostates. Each data point shows 
results for one macrostate. Two notable features are ap- 
parent. First, for isostatic states, the widths of the distri- 
butions AR^/ (Rl) increase with increasing; N. Second, 
the most compact packings for A^ = 10 and 11 are not the 



TABLE I: Numbers of macrostates M, macrostates with 
liquid-like (Ci) symmetry Mci , macrostates with Barlow or- 
dering M Barlow, and fractions of microstates with Ci sym- 
metry and Barlow ordering {fci and fsariow, respectively) 
^67j . Results for M for N < 10 are the same as reported in 
Ref. [25]. Note that some Barlow-ordered nuclei can have Ci 
symmetry, i.e. fsari + /ci > 1- *• A4(ll,27) excludes the 
"bridge" packings described in Appendix [A] 
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maximally contacting packings. For example, 25 of the 
259 A^ = 10, A^c = 24 packings have a smaller R^ than 
the most compact of the 3 A^ = 10, A^c = 25 packings, 
while for A/" = 11, the most compact Nc = 28 packing and 
the 66 most compact Nc = 27 packings have Rg below 
that of the Nc = 29 minimal energy packing. 
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FIG. 9: (Color online) Distributions of (a) and (b) shape 
anisotropy As- Results from left to right indicate increasing 
N and Nc. Roman numerals (I)- (IV) indicate the extremal 
packings shown in Fig. [12] 

Figure ^h) shows the shape anisotropy As for the 
same set of macrostates, i. e. As{N^Nc) = \J\maxl^min 
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where Xmax are the maximum and minimum eigenvalues 
of the moment of inertia tensor R^: 



10 



r1 



1 



N 



i=l 



(15) 



Here ei is the unit vector along the /-axis, where / = 
X, 7/, and z. Maximally symmetric (sphere-like) packings 
have As 1. Isostatic packings show a broad range of 
anisotropy that increases with increasing N . Anisotropy 
does not systematically increase with the degree of hyper- 
staticity. However, it is clear that the range of anisotropy 
decreases. Both the most and least symmetric packings 
are isostatic. 

To check whether the results in Fig. [9] are represen- 
tative of the full ensemble of packings, we examine the 
probability distributions P{R^q) P{^s)i where 



P{Rl){N,N,) = a 



M{N,Nc 

k=i 



(16) 



and 



P{As){N,N,) = n- 



MiN,N,) 
k=l 



s,ki 



(17) 



where R'^ ^s,fe, and ujk are the squared radius of gy- 
ration, anisotropy, and permutational entropy (Eq. [5| of 
the k^^ macrostate. 

Figure Il0| shows P{Rl) and P{As) for TV 



27. 



11, N, = 

SligTrny narrower distributions are obtained for 
smaller but are qualitatively similar. Results indicate 
that the most compact and the most symmetric nuclei 
have low entropy and are consistent with earlier stud- 
ies for smaller N [TT1|22]. P{As) is particularly broad. 
These distributions indicate that the "typical" nuclei is 
neither spherical nor characterized by a single value of 

Rg. 

Although widths of the distributions of shapes, sizes, 
symmetries, and entropies decrease for hyperstatic nuclei. 



they remain broad. The top panel of Fig. 11 shows the 
20 macrostates for TV = 11, Nc = 28; R^ and Ag are the 
same as in Fig. [9] The bottom panel shows their relative 
permutational entropies ujk/^ (Eqs.[5][6|. 

The structural diversity of packings reported above il- 
lustrates a key feature that should be included in theo- 
retical treatments of nucleation in hard- and sticky-hard 
sphere systems. Figure [TT] shows why one would not ex- 
pect classical nucleation theory to work for sticky hard 
sphere packings in this A^-regime and it is neccessary to 
consider nuclei possessing arbitrary geometry. There has 
been great interest in recent years in finding the dens- 
est finite sphere packings O [8], i.e. the A'-sphere nu- 
clei that minimize volume V. Most studies {e.g. Refs. 
[H El |43| ) search for nuclei optimizing either density or 
energy. For hard- and sticky-hard spheres, however, it is 
far from clear which quantity one should optimize. We 
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FIG. 10: (Color online) Probability distributions of (a) Rg 
and (b) shape anisotropy As for A = 11, Nc = 27 packings. 



have shown in detail that the relation between the num- 
ber of contacts and density is nontrivial, and that pack- 
ings optimizing these two features are in general different 
from each other. While this competition will break down 
at large A', i.e. the FCC crystal is simultaneously the 
densest packing and a maximally contacting packing for 
N ^ oo [1 , and the A"* at which the crossover occurs is 
unknown, we have shown that A'* > 11 [65 . 



Figure 12 further illustrates the competition between 
energy minimization and density maximization, and the 
importance of asphericity therein, by contrasting "ex- 
tremal" packings for A" = 11. The unique sticky hard 
sphere ground state (panel (I)) has Nc = 29 and HCP 
order [22] but is neither the densest, most compact, nor 
most symmetric. Panel (H) shows the packing that fits 
within the smallest spherical volume. It would be an 
ideal nucleus for hard-sphere crystallization within the 
framework of classical nucleation theory, but is in fact a 
"bad" (off-pathway) nucleus since it is a partial icosohe- 
dron lacking Barlow order. Panel (III) shows the densest 
packing (in the sense that it fits within the smallest con- 
vex shrink-wrapped volume Vgw)^ that has FCC order. It 
is also the most symmetric packing, i.e. the nucleus that 
minimizes Ag . Finally, panel (IV) shows the packing that 
is simultaneously the least compact and least spherically 
symmetric; it is also a "bad" nucleus lacking Barlow or- 
der. Note that the packings shown in panels (II-IV) are 
all second excited states for sticky spheres, energetically 
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FIG. 11: (Color online) Top panel: (a-t) The 20 distinguish- 
able macrostates for N — 11^ Nc = 28. Bottom panel: En- 
tropic fractions cj^/Q for these macrostates. Relative values 
of cjfc are inversely proportional the symmetry numbers c of 
the associated point groups [11 , i.e. the Ci macrostates (a-k) 
have c = 1 while the most symmetric macrostates have c = 8. 




FIG. 12: (Color online) The (I) ground, (II) most spherically 
compact, (III) densest, and most symmetric, and (IV) least 
compact and symmetric packings for N = 11. Packings (II- 
IV) are all isostatic (Nc = 27) and therefore second excited 
states. Sphere radii in panels (I-II) have been reduced for 
clarity. 



degenerate, and (except for (IV)) mechanically stable. 



B. Nuclei with structural motifs incompatible with 
bulk crystallization 

Stacking faults and five-fold symmetric structures are 
"defects" incompatible with bulk crystallization at ^ = 
n/y/lS. Their presence is well-known to impede hard- 
sphere crystallization [T2 Hl4| [20 ] . In this subsection, we 
quantify the propensity of nuclei to contain these and 
related structural motifs. 

The simplest stack- faulted motif is the M = 6, Mc = 
12 capped trigonal bipyramid structure shown in Fig. 
[TJb). Table [n] shows the number of macrostates Aictb 
and fraction of microstates fctb that include this motif. 
Note that Mctb and fctb are lower bounds for the num- 
bers and fractions of stack-faulted structures since other 
stack-faulted motifs exist. Nonetheless, the propensity 
for stack- faulting is surprisingly high given the small size 
of the nuclei - above 50% for all > 8 packings with 
Nc < N^^^. Stack faults appear in hyperstatic nuclei 
at N = 11, which is consistent with the fact that these 
nuclei do not maximize contacts. 

Two five-fold-like structural motifs often found in small 
nuclei are seven-sphere minimal energy packings {Mc = 
15) and are shown in Fig. 13 Panel (a) shows a five- 



fold-symmetric partial icosohedron. Panel (b) shows an 
nearly 5-fold symmetric structure which differs from (a) 
in that the green (lightest shaded) dimers contact and 
the 5-ring is open rather than closed with a separation 
rij = 1.19 [22 . Both of these structures are incompatible 
with close-packed crystal structure. 

FIG. 13: (Color online) Some structural motifs hostile to bulk 
crystallization: (a) The M = 7, Mc = 15 partial icosohedral 
and (b) M = 7, Mc = 15 open- loop near-5-fold symmetric 
structures. 



Table [ll| shows the fraction of macrostates for 7 < 
N < 11 containing motifs that are incompatible with 
long-range crystalline order. For isostatic packings, the 
5-fold symmetric subclusters in Fig. [13] are found in 
many macrostates and about 60% of microstates for 
9 < N < 11. In contract, these 5-fold symmetric sub- 
clusters are found in only one of the 24 nonisomorphic 
hyperstatic packings (and in < 1% of microstates). The 
near-5-fold symmetric structure shown in Fig. 13 ^b) cor- 
responds to an elementary twin defect. The high frac- 
tion of these (/sb) may explain why five-fold symmet- 
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TABLE II: Structural motifs that are incompatible with bulk crystallization as a function of N and Nc. Left columns: Propensity 
of nuclei to contain minimal stacking faults. Note that fctb is strictly > f^b since the stack-faulted structure shown in Fig. 
[T]^b) is a subset of the near- five- fold symmetric structure shown in Fig. 13 'b). Middle columns: Number of macrostates and 



percent of microstates containin g pa rtial-iscohedral five-fold symmetric (Mba and /sa) and open-loop {Ai^b and f^b) near-5-fold 



symmetric substructures in Fig. 13 Right columns: Propensity of nuclei to contain second nearest neighbors with < \/2; 
values of r^nd^ numbers of macrostates and fractions of microstates with r^^d ^ V^- Values of r^^^ 

for AT = 5, 

6, 7 and 8 are the same as in Ref. [22 . Note that there are several "near-miss" N — 11^ Nc — 27 macrostates with r^^^ only 
slightly above 1. Some of these possess soft modes, but we have verified that these cannot form 28th contacts. 



N 




Mctb 


fctb 




f5a 


M5b 


fdb 


^2nd 


M<V2 


f<V2 


7 


15 


3 


.286 


1 


.102 


1 


.153 


\/2 - 2/^5 


2 


.255 


8 


18 


7 


.567 


1 


.153 


3 


.166 


4^6/9 


5 


.331 


9 


21 


30 


.700 


7 


.098 


22 


.512 


\/2 - 2/^5 


33 


.691 


10 


24 


165 


.643 


32 


.135 


110 


.415 


1.03296 


185 


.721 


10 


25 




















V2 








11 


27 


1126 


.723 


220 


.130 


726 


.467 


1.00489 


1332 


.835 


11 


28 


8 


.511 








1 


.0116 


4^/9 


2 


.035 


11 


29 




















V2 









ric twinned crystallites are commonly observed in sticky- 
and hard-sphere systems [121 HO] • 

Another metric for nuclei incompatible with bulk crys- 
tallization at ^ = tt/vT^ is the minimum 2nd- nearest 
neighbor distance t^2' 

FCC-, HCP-, and Barlow- 
ordered crystallites have r^^2 — V^- Therefore nuclei 
with N > 6 and < cannot have Barlow order 

[66] . Table [n] shows the numbers of macrostates 
and fractions of microstates f^^ with r'^^ < V^- The 
fact that /<y2 + f Barlow < 1 for all A^ and Nc does not 
indicate any inconsistency since stack-faulted structures 
tend to be associated with neighbor distances r^j > ^/2, 
In the above subsections, we have examined several 
measures of microstructural order, f Barlow is a measure 
of "good" nuclei that are consistent with long-range crys- 
talline order (LRCO) while 1 - f Barlow, fba + /sb, fcth, 

and /<y2 independent measures of "bad" nu- 

clei that are inconsistent with LRCO [67 . The latter 
four all show the same trends; they increase with A^ for 
isostatic packings (to very high fractions) and decrease 
with increasing hyperstaticity. High energy barriers are 
expected between "bad" and Barlow-ordered nuclei since 
many bonds must be rearranged to change from one or- 
dering to the other. These results provide an explanation 
for the propensity of sticky hard sphere systems to jam 
and glass- form in both simulations and experiments. 



C. Applicability to other potentials and methods 

Our results for sticky hard-sphere packings are rele- 
vant for analyses of clusters formed by systems that in- 
teract via other potentials U{r) with hard-core like repul- 
sions and short range attractions. The sticky hard sphere 
model has been shown to provide a perturbative "refer- 
ence state" ^8] for such potentials (the sphere diameter 
D may be replaced by the minimum of a general inter- 



particle potential). In other words, packings of sticky 
spheres become increasingly similar to equilibrium clus- 
ters of particles interacting via a potential U{r) as the 
range of interaction Vc {i- e. U{r) = for r > Vc) ap- 
proaches D from above, and are rigorously identical in 
the limit of hard core repulsions and Vc/D 1. 

The minimum 2nd-nearest neighbor distance r^*^ (Ta- 
ble Il| is a particularly useful metric for evaluating the 
sticky hard sphere model's suitability for determining 
minimal energy clusters of other potentials since addi- 
tional local free energy minima begin to appear when 
rc > 1 + {r^fri - [43 . Our complete set of sticky 
hard sphere packings form an arguably complete set of 
initial guesses for identifying all isostatic clusters of up to 
9 particles when the interaction range {vc/D — 1) ^ .025, 
and for up to 10 particles when the interaction range 
{tc/D — 1) ^ .015. The hyperstatic packing sets should 
be suitable initial guesses for larger Vc] for example, 
Vc/D ^ 1 + (a/2 — l)/2 for our minimal energy (25— and 
29— contact) packings for A^ = 10 and 11, respectively. 

For larger Tc, "initial guesses" for strain- free iiiiiii- 
mal energy clusters can be obtained by selecting a subset 
of the sticky sphere packings satisfying 1 + (r^^^ ~ l)/2 ^ 
Vc [69 . Given a complete set of initial guesses for the 
set of nonisomorphic clusters and their permutational 
entropies {i.e. {r} and {uj} for each of the M{N,Nc) 
nuclei), sophisticated energy-landscape and transition- 
state analyses useful in cluster physics may be performed 
[HI [43] . Potentials for which such a procedure should be 
applicable include the "narrow" square well, the short- 
range limit of the Asakura-Oosawa and Morse potentials 
[11, 43 , and the hard-core attractive Yukawa potential 
in the strong screening limit. Such potentials describe 
a broad range of physical systems ranging from colloids 
interacting via depletion- mediated attractions [TTJ [20] to 
buckyballs [70 . Further, the ground and mechanically 
stable excited states of systems interacting via Eq.[l]have 
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been shown to describe the structure of real colloidal 
crystallites in dilute solution at ksT 46 [TTJ [TTJ [73] . 
These results provide justification for the use of Eq. [T]in 
our exact enumeration studies. 

Additionally, our results should be useful in numerical 
implementations of modern theories for aspherical nu- 
clei based on cluster expansions (Mayer f-bond diagrams 
j72]). Such models are commonly used in liquid-state 
theory [iTlEHlllI], and are in principle exactly soluble, 
but have to date suffered from incomplete sets of Mayer 
diagrams describing differently structured aggregrates as 
well as "dangerous" [27 singular cluster integrals that 
impede implementation of such theories. The contact 
graphs {i.e. valid adjacency matrices) reported here [64] 
correspond to a complete set of Mayer f-bond diagrams 
for ground state aggregrates of M < 11 particles as well 
as stable first excited states (for M = 10) and both first 
and second excited states for M = 11, while the implicit 
contact graphs (Appendix [C| correspond to cluster inte- 
grals that are singular. 



ergetic cost. 

In relating our studies to crystal nucleation we assume 
that mechanically stable nuclei play a key role. This 
claim is clearly well supported in the dilute regime where 
solvent effects are minimal [TTl |75]. Additionally, our 
contention that understanding the statistical-geometrical 
properties of small nuclei (and in particular, the promi- 
nence of structural motifs that are incompatible with 
the bulk crystal) is relevant to glass-formation and 
jamming is consistent with recent experimental work 
by Royall et.al [62], which indicate that such local 
motifs lead to kinetic arrest in colloidal suspensions 
possessing hard-core repulsive and short-range attrac- 
tive interactions. In the semidilute or concentrated 
regimes in which solvent effects are stronger, while our 
study cannot capture all the complexities of nucleation 
from the bulk, quantitative comparison of the nuclear 
structures reported here to those reported in studies of 
crystallizing liquids [16] [171 El] is an interesting topic 
for future work. 



IV. DISCUSSION 

An advantage of exact enumeration is that it identi- 
fies all nuclei that can form, as opposed to only those 
that do form for a specific preparation protocol. How- 
ever, it treats A/'-particle nuclei in isolation and neglects 
solvent effects. Sticky hard-sphere nuclei in a solution of 
other sticky hard spheres would be "continuously bom- 
barded by and grow by absorbing smaller clusters" [T9] . 
Such collisions can infiuence the pathways by which small 
nuclei form larger crystallites. For example, an "off- 
pathway" (non-Barlow) nucleus might be excited by a 
collision and reform into a larger "on-pathway" Barlow 
nucleus. Additionally, the detailed structure of nuclei in 
such a solvent would be altered both by finite tempera- 
ture {e.g. vibrational entropy pJJlTSJ) and the crystallite- 
fiuid interfacial free energy [29]. 

Many Monte Carlo studies have examined crystalliza- 
tion in bulk hard-sphere systems p!2HT5] . Schilling et.al. 
[15] performed Monte Carlo simulations of crystallization 
in dense {(j) = 0.54) hard-sphere liquids and argued that 
crystallization occurs through a two-step process wherein 

(1) dense "amorphous" clusters form and act as (2) "pre- 
cursors" for nucleation of larger close-packed crystallites. 
They identified the growth of Barlow order during stage 

(2) using the bond-orientational order parameter q^, but 
did not examine the detailed structure of the amorphous 
clusters {i.e. they did not examine stacking faults or 
5-fold symmetric structures) [74 . Our results are not 
fundamentally inconsistent with theirs; the non-Barlow- 
ordered nuclei we have identified above could correspond 
to their "amorphous" clusters. We also note that dif- 
ferences between sticky-hard and purely repulsive hard 
spheres will significantly alter the physics of the (1)^(2) 
process since non-Barlow clusters must break bonds to 
rearrange into Barlow clusters, with a corresponding en- 



V. CONCLUSIONS 

In this manuscript, we described the structural prop- 
erties of a complete set of isostatic and hyperstatic pack- 
ings for hard spheres obtained via exact enumeration. 
For sticky hard spheres with contact attractions, we also 
analyzed mechanical stability. Our key findings included 
exponential growth in the number of nonisomorphic iso- 
static packings and nontrivial variation of the size and 
symmetry of packings with increasing hyperstaticity. We 
also calculated the absolute and relative entropies of all 
packings and their propensity to include various struc- 
tural motifs that are either compatible or incompatible 
with bulk crystallization at (/) = tt/a/TS- Isotatic nuclei 
form an increasingly liquid-like ensemble as N increases. 
For example, the fraction / of isostatic nuclei possessing 
Barlow-order decreases rapidly with N to only about 5% 
for TV = 11, and the remaining 95% contain defects such 
as stacking faults and 5-fold-symmetric substructures. 
While / increases with hyperstaticity H = Nc — Njso^ f 
is only about 50% for A/" = 11, H = 1 nuclei. Although 
we terminated our enumeration studies at = 11 due to 
the limits of current computational resources, the trends 
reported here should [H |12] continue to hold for higher 
N. 

Additionally, we have shown that considering nuclei 
with N ^ 10 captures a complex regime [65j where clas- 
sical nucleation theory performs particularly poorly [18] . 
In this regime, maximizing density and maximizing Nc 
compete, the distributions of nuclear size and symmetry 
are broad, and many nuclei are highly aspherical. Since 
colloids with hard-core-like repulsions and short-range at- 
tractions form stable nuclei in this size regime [TTJ [76] , 
our results present challenges for traditional theoretical 
approaches to nucleation in sticky hard sphere and re- 
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lated systems. Most analytic and seminumerical treat- 
ments {e.g. phase field theory and the classical density 
functional of Cahn and Milliard [78 ) either assume or- 
dering consistent with the bulk crystalline phase, or al- 
low for ordering different than that of the bulk crystal 
but assume that nuclei are spherical. Our results sug- 
gest that such restrictions in traditional methods prevent 
them from capturing the potential complexity of small- 
N crystallite nucleation. Novel theoretical treatments of 
nucleation should consider nuclei of both arbitrary order 
and arbitrary geometry. 

Finally, while there have been many recent detailed 
studies of crystallization in hard-sphere systems, there 
have been relatively few [71 ISTJ [75] theoretical studies of 
the dynamics of sticky hard sphere crystallization. The 
higher relative entropies of cluster formation for less- 
ordered nuclei should strongly affect nonequilibrium be- 
havior. For example, the large fractions of small nuclei 
with Ci {i.e. liquid-like) symmetry, fivefold symmetry, 
stacking faults, and other types of non-Barlow ordering 
constitute an effective "entropic" barrier to nucleation 
and growth of large ordered crystalline domains that 
should play a key role in controlling the critical quench 
rate above which these systems glass-form. For sticky 
spheres, mechanical stability of the non-Barlow nuclei 
presents an additional energetic barrier to ordered crys- 
tallite growth. It would be interesting to compare the 
ensembles of nuclei produced in nonequilibrium studies 
of sticky sphere aggregation to "ideally prepared" ensem- 
bles (in which all possible aggregates are obtained) like 
those reported in this paper. 
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Appendix A: Mechanical stability: Floppy packings 
and bridge structures for = 11 

"Floppy" packings possessing soft modes are of spe- 
cial interest for nucleation studies at finite T since they 
possess higher vibrational entropy [11] . Table |III| shows 
values of the number of macrostates and fractions of 
macrostates possessing "nontrivial" soft modes (i.e. soft 
modes in packings that satisfy the two necessary condi- 
tions (i-ii) for mechanical stability). These soft modes 



correspond to small collective monomer motions that do 
not break contacts [38 . Interestingly, the number of non- 
trivially floppy macrostates Ai floppy increases faster than 
A4 with increasing N over the range 9 < N < 11. Both 
this result and the exponential increase in M (Table |l| 
are related to the emergence of hyperstatic packings with 
Ar^^^^(lO) = 3^" - 5 and N^''^{11) = 3N - 4. The 
presence of these hyperstatic states both makes it eas- 
ier to form isostatic packings (by adding low-coordinated 
spheres to a hyperstatic packing) , and increases the like- 
lihood that such packings will be floppy. However, the 
latter effect is small for the range of N considered. 



TABLE III: Number of macrostates and fraction of mi- 
crostates possessing nontrivial soft modes. No N < 9 packings 
possess soft modes, and none of the < 11 packings possess 
more than one soft mode. 



N 


Nc 


M floppy 


/floppy 


9 


21 


1 


.00427 


10 


24 


4 


.0194 


10 


25 








11 


27 


31 


.0136 


11 


28 


1 


.0116 


11 


29 










FIG. 14: Color online) Structure of (a) trivially and (b) hy- 
perstatic floppy nuclei for = 11. Panel (a) shows a typical 
bridge packing; a bridge monomer possessing only two con- 
tacts is shown in green (the lighter shad e) . Panel (b) : This 
AT = 11, A^c = 28 packing ((q) in Fig. [iT]) is the smallest 
hyperstatic sticky sphere packing possessing a soft mode. Ar- 
rows indicate particle displacements proportional to a non- 
trivial zero eigenvector of the dynamical matrix. 

Since 10 is the smallest N at which hyperstatic pack- 
ings can form, it follows (but has not heretofore been 
shown) that 11 is the smallest N at which "bridge" struc- 
tures can form. These structures are "trivially" floppy 
because they include monomers that possess only two 
contacts and fail to satisfy condition (i). We find 25 
graph- nonisomorphic A/" = 11, A/'c = 27 bridge packings. 
A typical example is shown in Fig. [l4ja); bridge spheres 
(shown in green) have a single floppy mode associated 
with the free configurational degree of freedom (motion 
along a circle of R = V3/2 centered on the line con- 
necting the two contacted spheres). Since we focused on 
mechanically stable or nontrivially floppy nuclei, these 
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25 packings are not included in our structural analyses 



in Section [nil However, it is important to include bridge 
packings in exact enumeration studies that consider fi- 
nite temperature. These structures possess high config- 
urational entropy, and for sticky sphere systems, should 
dominate equilibrium populations of = 11, A^c = 27 
nuclei when ksT is not small compared to the contact 
energy e. 

Another interesting feature of = 11 packings is that 
11 is the smallest N at which a hyperstatic sticky sphere 
cluster with a floppy mode can form. The packing and 
associated floppy mode are illustrated in Fig.[l4]^b). This 
packing has topology similar to a subset of a BCC lattice; 
3 adjacent squares of bonds surround a linear trimer. The 
floppy mode is a torsional motion about this trimer, and 
is associated with the "squares" of bonds forming the 
outer part of the packing. 



Appendix B: Verification of A4 for N = 11 




all cases, including earlier tests on = 10 packings (not 
shown), convergence is found for ut ^ 100 and increasing 
ut by another order of magnitude produces no additional 
solutions. Faster convergence is found for smaller N. 
Appendix C: Implicit Contact Graphs 

Equation [S] is a set of equations and inequalities suffi- 
cient to obtain the structure of A/'-sphere packings with 
at least Nc contacts and no overlaps. The geometric re- 
jection rules enforce only the ">" portion of the Tij > 1 
condition for {i^j} pairs with Aij = 0. For N > 10, 
there exist A with Nc contacts whose Euclidean solution 
{r} is a packing with Nc -\- 1 contacts. Such "implicit 
contact" graphs violate the spirit of our enumeration 
method. Therefore, all A containing M x M implicit- 
contact subgraphs are rejected at the M-rule level. 

We find 29 nonisomorphic M = 10, Mc = 24 implicit- 
contact graphs. All reduce to one of the three M = 10, 
Mc = 25 minimal energy packings [22 when we solve for 
{r}. Note that we overcounted A^(io,24) in Ref. [37] 
by including 20 implicit-contact macrostates. Similarly, 
we find 13 M = 11, Mc = 27 graphs, not containing 
isomorphic subgraphs of the abovementioned set of 10/24 
graphs, that imply structures with 28 or more contacts 



Fig. 16 



10 20 



100 200 



illustrates one such structure: panel (a) 
shows an implicit contact adjacency matrix with A^ = 
11 and Nc = 27; contacts present in A are shown in 
black. The corresponding packing (panel (b)) possesses 
29 contacts (as determined by solving for {r}) and is 
identical to the N = 11, Nc = 29 packing shown in Fig. 
[T2|^a). The implicit contacts are shown in red in panel (a). 
In panel (b), the implicit contacts are shown as dashed 
lines, and spheres possessing implicit contacts (particles 
4, 6, 8, and 10) are shaded blue (dark). 



FIG. 15: (Color online) Convergence of structure-sol ver fo r 
N = 11 packings with increasing nr (see Section 
M final = 20 for A^e = 28 (dotted line) and 1658 for Nc 
(solid line; this includes the 13 implicit contact packings dis- 



IID). 
27 



(a) 



cussed in Appendix [c]). Here ni — 40 and TOL 
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For A^ > 9 our set of rejection rules is incomplete and 
we rely upon the structure solver to determine whether 
packings are valid. Here we provide evidence validating 
this procedure by showing convergence with increasing 
Ut of the fraction of adjacency matrices identified as valid 
(see Fig. [I]) Fi sure 15 shows results for a set of 43 non- 
isomorphic A^ = 11, A'c = 28 packings and a set of 4534 
nonisomorphic N = 11, Nc = 27 packings that passed a 
set of M < 10 rejection rules from a non-final version of 
our code. For fixed n^, the structure solver finds M solved 
valid packings. Msoived converges to M final in the limit 
of large ut and we report Ai = Ai final 
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FIG. 16: (Color online) Panel (a): An implicit contact adja- 
cency matrix with M — 11, 27 explicit contacts (black I's) 
and two implicit contacts arising from the solution of Equa- 
tion [s] (red [lighter shaded] I's). Panel (b): Visualizing the 
structure reveals a A^c = 29 packing; implicit contacts are 
shown as dashed green lines. 
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